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This paper will devote to construct a family of fifth degree cubature formulae 
for n-cube with symmetric measure and n-dimensional spherically symmetrical 
region. The formula for n-cube contains at most + 5n + 3 points and for 
n-dimensional spherically symmetrical region contains only + 3n + 3 points. 
Moreover, the numbers can be reduced to n'^ + 3n + 1 and n^ + n + lifn = 7 
respectively, the later of which is minimal. 
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1 Introduction 

In this article, the fifth degree formula with respect to the following integral 



will be considered where the region C M" and p(X) > for X G We study two 
cases: one of which requires Vt = [—1, 1]" and p(X) = pi(xi) . . .p„(x„) where pj(xj) is 
symmetrical; and the other of which requires that =^ is a spherically symmetric integral 
(see the definition of section 3 or refer to [9]). 

Among all the cubature formulae, those with few points received most attention. 
Without any special assumption, a general lower bound was given as follows (see [2] and 
the reference therein): for a cubature formula of degree 2k or 2k + 1, the number of the 
points > dimPj^, where denotes the polynomial space of dimension n with degree 
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Abstract 




(1) 



no more than k. For even degrees, no greater lower bound was given until now, see 
Cools' survey paper [2]. However, for odd degrees 2k + 1, the lower bound for centrally 
symmetrical integral can be improved as follows (see Moller [TT]): 



N > 2dimP^ - 



if s is odd, 

1 if s is even. 



(2) 



or explicitly. 




N > < 



V 




+ 5^(1-2—) 



s=l 



s + k-1 



s 



if k is odd. 



if k is even. 



(3) 



Particularly, for the fifth degree case 



> + n + 1. 



However, there is a big gap between the theoretical lower bound and the numbers of 
points in the known cubature formulae, especially for the higher dimensional case (see 
Cools' survey papers O 01 16] ) . So far, only for special dimension, the minimal cuba- 
ture rules can be obtained. Noskov and Schmid [13] derived a necessary and sufficient 
condition of existence of the minimal cubature rules with respect to the integrals over 
the ball. 

For the product region, a good choice to construct cubature formula is to employ 
Smolyak's method [H [TJ [HI [15] . However, the number of the points is much greater than 
the theoretical lower bound. Very recently, starting from Smolyak method, Hinrichs 
and Novak [8] presented a kind of formula with + 7n + 1 points for fifth degree 
case and (n^ + 21n^ + 20n + 3)/3 points for seventh degree case. Before them, except 
Lu's formula the best upper bounds were of the form "r; 2n^" and "~ 4n^/3" 
respectively, where "~" denotes the strong equivalence of sequences. Lu's formula uses 
n"^ + 3n + 3 points just for entire n-dimensional space with Gaussian weight function, 
which is the minimal point formula for the general dimension among all the known 
formulae until now. 

In this paper we shall present a method to construct fifth degree formulae with 
at most + 5n + 3 points for the integrals with the product form and with at most 
n"^ + 3n + 3 points for the spherically symmetrical integrals, which employ fewer points 
than Hinrichs and Novak's formulae [8] and are more close to the theoretical lower 
bound. This paper proceeds as follows. Section 2 will present a method to construct 
fifth formula for n-cube with symmetrical product measure. And section 3 will devote 
to the construction of cubature formulae for spherically symmetrical integrals. Finally, 
section 4 will give some numerical integration rules by our method. 
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2 Fifth degree formula for symmetrical product mea- 
sure 

We begin with a well-known formula with respect to the surface of sphere: 
t/„ = {X e M" : X? + + . . . + = 1}. 
Due to Stroud [16], the monomial integrals for p{xi, X2, ■ ■ ■ , Xn) = 1 are 



V 



Ur 



"r(n/2)' 

x-x- x-dx = 2r[K + i)/2]...r[K + i)/2] 

^ 2 ••• n r[(n + ai + a2 + ... + a„)/2] ' 



provided all the are even; otherwise the integral is zero. 

Mysovskikh [12] derives a cubature formula of degree 5 (n > 4) as follows: 

„ n+l n(n+l)/2 

ifi(/)=/ /(X)rfX^g(/) = A^[/(a(^))+/(-a(^-))]+5 5^ [/(b(^))+/(-b(^) 

i=i j=i 

(4) 

aM = 4^\...,a(:)), r = 1,2,..., n + l, (5) 



where 



_n+l ^ ^ ^ 



n(ji-i+2)(n-i+l) ' 



«f^ = ^/^^^±i^i^, ^ = r, (6) 

n(n—r+2) ' 

0, i > r. 



and 



{b(^) = {b?, hf, . . . , 6(f))} = { + a^'^) : A: < /, / = 1, 2, . . . , n + l}. (7) 

The corresponding coefficients are 

A = n{7 -n)V/[2{n + lf{n + 2)], 
B = 2(n- l)V/[n(n + l)2(n + 2)]. 

We firstly construct a cubature formula of degree five for the integral with product 
form: 



^(/)= / Pl(Xi)p2(x2)...Pn(x„)/(X)dX, 



where 



J pi{xi)dxi = 1 and j pi{ 



(xi)xf +Mx, = 0. 



Define a diagonal matrix T — diag{aii, 022, • • • , (inn\-i where Oj 



a positive parameter 7. Let a^-'^ = Ta^-'^ b*^*^^ = Tb*^'^^ A = 7^4, B = 7^ and 

n+l n(n+l)/2 

j=i 3=1 



Then 



Q{xi) = 



r=l 



n+l 



r=l 

n(n+l)/2 



r=l ^=1 



■ ^^i{XiXj) 



7i^i(x?xi) 
f^PML.^(.|), for 
[0, for i^j] 

^^•V7-^(-D, for i = 



0, 



for i ^ j; 



with 



(9) 



n+l 



Q(xiXkXi) = ^ [a^af ^ • akka^^^ ■ aua'f^ - aua'f^ ■ akta^^'' ■ maf'^] 



(0 „ Jr) 



>(0 „ Jrh 



r=l 

n(n+l)/2 



(r) 



r-=l 



n+l 



Q{xiXjXkXi) = ^ E [^Onol''^ • ajja'f^ ■ akkat^ ■ aiiat^] 



r=l 

n(n+l)/2 



n(n+l)/2 



r=l 



n+l 



r=l r=l 



J!f{xfxi), for i = j and k = I; 

i^i^)f-^V for ^ = J = k = h 
0, others 

■^{^l^Di fo^ ^ = j and /c = /; 
3 {^{x'i)f , for i=j = k = l; 
0, others 

and similarly = 0. Define ^ = ^ — Q and compute the moments as 

follows: 

^(l) = l-7-^, 
^{Xi) = 0, 

{0, for i ^ j, 

- 3 {^{x))Y , for i = j = k = I, 
0, others, 

{XiX jXl^XlXm) 0* 



=Sf {XiXjXf^Xi^ 



Remark 1. If R{f) is a cuhature formula of degree five with respect to S£ , then =Sf ~ 
Q + R is a cuhature formula of degree five with respect to . 

Now it is the time to state our main result. 

Theorem 1. There exists a cuhature formula Cn for =Sf of degree five which uses at 
most ra^ + 5n + 3 points. 

Proof. Suppose that an integration formula R of degree five for ^ is given by 
v(*'^) = Vi^jGi Wij, z = 1, 2, . . . , n, J = 1, 2 

(0,0. ..,0) C ^ ' 

where e^ = (0, . . . , 0, 1, 0, . . . , 0) is a vector of dimension n and the ith element one is 
the only nonzero element. To enforce polynomial exactness of degree five, it suffices to 
require (fTO|) to satisfy the following equations: 

2 

w^, {v,^,r = A^T) , m = 1, 2, 3, 4, 5, z = 1, 2, . . . , n, 

(11) 

n 2 ^ ' 

5^5^^,, + c = ^(i). 

i=l j=l 
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By a proper parameter 7, we can make ^{xj)-^{xf) > 0. Thus we get the required 
cubature formula R with 



n 

1=1 

The corresponding cubature formula can be written as 

n+l K(n+l)/2 
n 2 

+EE^^.^-^^'''^+^/(0'0'---'0)- 
i=i j=i 



(12) 



(14) 



The total number of the points is at most [r? + 3n + 2) + (2n + 1) = + 5n + 3. □ 
Remark 2. To ensure all the points inside the integration region, it is required that 

au<l, and < < 1. (13) 

If^ix^j) < 3=^(x2)2, then 

7>2.-t^(n/2 + 2)(^(.D)^ V7>^^g|g5), 

(if (xf ) - ^{xt) + 3^(x,^)^) V2r (n/2 + 2) 
- JSf (xf ) (n + 2) W4 

which always has a solution. Otherwise 

(^(x^ 2r(f + 2) 

^i(x?xi) - 7rt(n + 2)' 

usually does not derive solution which means some points lie outside of the region. 
Particularly, if pi{xi) = (1 — x?)"^*(Q!i > —1), then 

6 

^{x',) - 3^(xJ)^ = "(a,/2 + 3)2(a,/2 + 5) < ° 
implies that the solution for 7 always exsits. 

Remark 3. 5'mce A vanished when n — 7, the total number of the points is + 3n+l 
ifn = 7. 
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Remark 4. // 

y2vrf(n + 2) 1(1 + 1) 

2v/r(f + 2) ' ^ 7r-/2(n + 2)' 

t/ien =Sf = anc? thus 

7-Q(/) + (l-7-^)/(0,0,...,0) 

is a cubature formula of degree three. Furthermore, it is easy to check that 

7■Qif) + J2l^^f{0,0,...,0) + il-J2l^-7■V)f{0,0,...,0[ 

1=1 * i=l 

is a cubature formula of degree five provided f is smooth enough, where 



4! 



3 Fifth degree formula for spherically symmetrical 
measure 

^ is a spherically symmetrical integral functional means the region is of the form 
{X G W^\p < ||X|| < q] and the weight function W{1C) is a function of ||X||, where 
< X, Y >= xiyi + . . . + XnUn and | |X| | = v^< X, X > (see [9]). The classical examples 
are the integrals over the ball and over the entire space with Gaussian weight. In 
this case, .^{xf) = =Sf(xj) and ^{x^xj) = ^{x\x'^) for any i,j,k,l. Taking an = 

^^^/f^L and the same mathod explained, we have 
Q{xi) = 

[0, for iy^j; 

Q{xiXjX}^ 0, 

{^{x}x\), for i = j and k = I] 
3 {^{xlxD) , for i = j = k = I; 
0, others. 

QiyXiXjXj^XiXfyi^ 



and 



^(l) = l-7-^, 
^{Xi) = 0, 
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0, for i ^ j, 

■^^"^^"^^■^ = ^(x?) - ^^§t^ . V7 ■ ^{xlxD, for ^ = J, 



2^r(f+2) 

'^(xf) -3(^(x2x2)), for i=j = k = l, 
0, others, 



where ^ = ^ - Q. 

It is a shghtly different from the product case, because we have the following 

Theorem 2. Suppose that 

^(/) = / /(X)p(X)dX 

is a spherically symmetric integral functional, then 

££{x\) - 3^(x?a;2) = 0. 

Proof. To calculate ^{xf) and ^{x\x\), we introduce the transformation which was 
used in [16], Page 33] 

Xi= r cos 6n-l cos ft„_2 • • • COS O2 COS 6'i 

X2 = r cos 6'„_i cos ^^„_2 • • • cos 6*2 sin 9i 

xs = r cos6'„_i cos6'„_2 . . . sin6'2 , ^, 

(15) 



x„_i = r cos 6'„_i sin 6'„_2 
= rsin6'„_i 

The Jacobian of transformation (ITS!) is 

J = r"-i(cos^„_i)"-2(cos^„„2)""'. • • (cos^3)'(cos02) 

and therefore 

^(xt) = I r+ / p{r)\r\"-\^dr [ ' (cos^i)M^i x Jq 
££{x\xl) = \ r + [ '')p(r)|rr"Wr /" ^ (cos^i)2(sinei)Mei x Jq 

Up J-q ) J-tt/2 



where 

/o= r^'(cos^2)M^2 r^\cosesfdes... r^\cosen-ir+Me.,-i. 

J-tt/2 J~-k/2 J--K/2 



Hence 

f7r/2 



^{xlxD /^/^2(cosei)2(sin^i)M^i 
and therefore ^{x\) — ?>^{x\x\) = 0. This completes the proof. □ 



Let 



=Sf (xf) , ' ■ \ H ■ £^{x\xl) = 



2v/r(i + 2) 

that is 

_ i^ixj)y 2T{n/2 + 2) 
^ ~ ^{xlxD ' 7r"/2(n + 2)2' 

In this case, we get a fifth degree formula with at most + 3n + 3 points 

n+l n(n+l)/2 

= i ^ [/(a(^-)) + /(-a(^-))] +B Yl [/(b^^'^) + /(-b^^'^ 
i=i 3=1 
+ (l-7-V^)/(0,0,...,0). 



(16) 



If n = 7, A = will yields a minimal cubature formula which only contains + n + 1 
points. 

4 Numerical Results 

Assume n > 4. 

• Symmetric product region. 

We consider the case of weight functions pi = p2 = • • • = Pn = 1/2. 
In this case, 

^{xt) = ^(xj) - 3(^(x2))2 = 1/5 - 1/3 = -2/15 < 0. 
To confirm all the points inside the region, it follows from ( 11^ that 

2r(n/2 + 2) 



7 



> 



If we take 7 = then 

V30(r2-1) 



5(r2- 1) 

9 



Wij = — 777^?^^ + —n — 777^, C = —n"^ — -^n^ — —n + 1. 



8 2 7 
— n 

108'^ ' 54'^ 108' 54'^ 27 54 

Due to (fT2!) . the corresponding cubature formula can be written as 



^ ^ /(X)dX^(^n3-|-n2-ln + l)/(0,0,...,0) 



2"./r_i,i]n^ ' '54 27 54 



+ ^ 108^ 54 108^^l^^l^ 5(n-l) ^'^J+^l^ 5(n-l) 
j=l 

• Spherically symmetric integral 

(1). Entire n-dimensional space with Gaussian weight function. 
Simple calculation yields 

r(n/2 + i) , ^- 

7r"/^(n + 2) 

Thus the fifth degree formula can be written as 

f _ 2_ _ 2 27r"/2 
/ e"'^i"-"'^"/(Xi, X2, . . . , Xn)dxi ...dxn^ — ^/(O' 0, . . . , 0) 
Jr" ^ + 2 

^ ,(: 1^).";!; g [/( ■ ^ /(- vw^- 



which is same with Lu's formula [TO] . 
(2). The n-dimensional unit ball. 
Simple calculation yields 



2(n + 4)r(n/2 + 2) , n + 2 
7 = 5 ■ , and an = \ -. 

(n + 2)V"/2 " Vri + 4 

Thus the fifth degree formula can be written as 

87r"/2 

f{Xi, X2, . . . , Xn)dXi ...dXn^ —. -^^^^^^T^ f ' 0, . . . , O) 

(7 - n) n (n + 4) 7r"/2 ^ /^^^ ^^^.^^ ^ . (i) 



(n + l)^(n + 2)''r(n/2)^' Vn + 4 " ' V n + 4 
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(« + 2)'r(n/2) ^ n + 4 " ' + V n + 4 " 'J' 

(3) . The n- dimensional Spherical Shell. 
The region S^^^^ is 

^^'^^^ = {(xi, X2, . . . , x„)|r < X? + + . . . + < 1} 

and weight function p(X) = 1. The monomial integrals can be found in ^16j. In this 

case, 

(n + 2)V"/2 Vn + 4 

Thus the fifth degree formula can be written as 

f 87r"/2(l - r'') 

fix,, x„ . . . , x^)dx, . . . dx„ - ^^ + 2)^r(n/2) ^(°' °' • ■ • ' °) 

I i7-n)n{n + 4) vr"/^(l - r») ^ r ./^. a^')) 

(n + l)'(n + 2)^r(n/2) \n + A 

n(n + l)^(ii + 2)'r(ii/2) j^i V" + 4 

+/(-^'l^.y|||.bC')]. 

The origin is outside of the region. 

(4) . Entire n- dimensional space with weight function exp{— xf + + • • . + x^). 
Simple calculation yields 

(n + l)r(n/2 + 1) , r 

Thus the fifth degree formula can be written as 

f / 2^ ^ 2 (2n + 3)7rV2"+2r(^4i) 

Jr" (n + 2)(n + 3) 



4(ri-l)27rV2"r(^^i^"^"+'^/' 
(n + 3)(n + l)(n + 2) 
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